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Calculations of high-energy processes involving the production of b quarks are typically per¬ 
formed in two different ways, the massive four-flavour scheme (4FS) and the massless five- 
flavour scheme (5FS). For processes where the combination of the 4FS and 5FS results into a 
matched calculation is technically difficult, it is possible to define a hybrid scheme known as 
the doped scheme, where above the 6-quark threshold the strong coupling runs with n/ = 5, 
as in the massless calculation, while the DGLAP splitting functions are those of the n/ = 4 
scheme. In this contribution we present NNPDF3.0 PDF sets in this doped scheme, com¬ 
pare them with the corresponding 4FS and 5FS sets, and discuss their relevance for LHC 
phenomenology. 


Motivation The accurate modeling of heavy quark mass effects is a crucial ingredient of LHC 
phenomenology. In processes involving 6-quarks, two calculational schemes are usually employed: 
the massless nj = 5 scheme (5FS), where the 6-quark mass is neglected and potentially large 
logarithms of the type InQ^/m^ are resummed into a 6-quark PDF b(x,Q‘^), and the massive 
rif = A scheme (4FS), where the non-zero value of mi, is kept in the calculation. Each scheme 
has its own advantages and weakness. 

A general framework for the combination of 4FS and 5FS results into a single matched 
calculation, FONLL, is available, and has been applied to several processes from heavy quark 
hadroproduction^ to DIS structure functions^ and Higgs production via bottom quark fusion^. 
The implementation of FONLL is however technically challenging for processes for which an¬ 
alytic expressions are not available (such as in semi-automated NLO codes) and therefore a 
hybrid scheme that captures the essence of this matching may be useful, even though it will be 
necessarily less accurate than the full matched calculation. 

As discussed in Ref. for many LHC processes the 4FS is perfectly adequate, since po¬ 
tentially large collinear logarithms in rrib are suppressed both by the smallness of the 6-PDF 
at large-x and by universal phase-space factors. However, a known drawback of this scheme 
is that with the 4FS evolution of the [3 function for the running of the strong coupling, above 
the 6-quark threshold Q = nih the value of as{Q) is not accurately reproduced, and therefore 
the strong coupling running cannot be simultaneously correct at relatively low scales, relevant 
for many of the datasets that enter global PDF fits, and at the higher scales relevant for the 
majority of LHC processes. This difference induces a suppression in the predicted cross-sections, 
which can be especially important for multi-leg processes at the Born level where the calculation 
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Figure 1 - Running of the strong coupling as{Q) in the 5FS and 4FS, normalized to the 5FS result. 


starts with a high power of 

To overcome these difficulties, a hybrid scheme for the treatment of heavy quark mass effects 
has been proposed^, named the doped scheme, where as{Q) runs with Uf = 5 flavors above rrih, 
thus having the standard value at Mz (consistent with the global PDG average), while DGLAP 
evolution is still performed in the Uf = 4 scheme, and in particular no 6-quark PDF is introduced, 
as in the massive scheme. This scheme allows the combination of the useful properties of the 
4FS, in particular accounting for bottom quark mass effects, while also using a value of as{Q) 
consistent with the global average for all values of Q. 

Therefore, the basic idea of the doped scheme is to use, above the 6-quark threshold, a 
four-flavor factorization scheme (therefore a 4FS for the DGLAP evolution) with a five-flavor 
renormalization scheme (thus, a 5FS for the running of as((5)). However, as explained in Ref.®, 
this choice is inconsistent, because then there would be a mismatch between terms in the splitting 
functions (specifically Pgg) which depend on the beta function, and those governing the evolution 
of the strong coupling. Recently, it was suggested in Refs. that this inconsistency can be cured 
by subtracting order by order from the partonic matrix elements the difference between the 4FS 
and 5FS running of as{Q), up to the order at which the hard cross-section is computed. Once 
this subtraction is implemented, one consistently restores the 4FS up to the finite order at which 
the matrix element has been computed. 


NNPDF3.0 sets in the doped scheme We now present results for a version of the NNPDF3.0 
NLO and NNLO global PDF sets® that are suitable to be used the doped scheme calculations, 
and compare them with their 4FS and 5FS counterparts. Following the standard method adopted 
in the NNPDF framework®, the doped sets are constructed starting from the NNPDF3.0 Uf = 5 
global fit at a scale equal to the 6-quark mass, = nil, then evolving upwards in using 
the DGLAP equations but now with nj = 5 in the /3-function for the running of asiQ"^) and 
with n/ = 4 in the splitting functions. The DGLAP upwards evolution is performed with the 
APFEL program^®, using the same truncated solution of the evolution equations as that used in 
the NNPDF3.0 fits. 

First of all, we show in Fig. 1 a comparison of the running of the strong coupling as{Q) in 
the 5F and 4F schemes. The different running in the 4F as compared to the 5F scheme induces 
variations on the value of as{Q) between 3% at Q ~ 10 GeV, and 6% at Q ~ 10 TeV, which 
in turns would lead to a shift 5a/a ~ nSag/ag for processes whose 4FS calculation starts as 
(T ~ a” at the Born level. 

In Fig. 2 we compare the gluon PDF, as function of Q, for two values of x, x = 10“^ and 
X = 0.1, at NLO and at NNLO, between the 5F (reference), 4F and doped schemes. As can be 
seen from the comparison, at large-x the doped scheme leads to intermediate results between 
the 4F and 5F schemes. On the other hand, at small-x, the doped scheme leads to a stronger 
rise of the gluon at small-x as compared to the 4F calculation, due to the larger value of as{Q) 























Figure 2 - Comparison of the NNPDF3.0 global fits in the 5F, 4F and doped heavy quark schemes, as a 
function of the factorization scale Q for two different values of a:, a; = 10“^ (upper plots) and a; = 0.1 
(lower plots). We show results normalized to the central value of the 5FS PDF set, for both NLO (left 
plots) and NNLO (right plots). 


at high scales, which drives the small-x DGLAP evolution. The results are qualitatively the 
same at both perturbative orders, NLO and NNLO. The same comparison, now for the down 
quark PDF d{x, Q^) is shown in Fig. 3. At large values of x, where effects of DGLAP evolution 
are moderate, the doped scheme coincides essentially with the 5F scheme, and the differences 
with the 4FS are a few percent at most. At small-x, the doped scheme also leads to a steeper 
rise of the light quark sea. 

Implications for LHC phenomenology A first study of the use of the doped scheme at the 
LHG for the production of vector bosons associated with a bb pair was presented in Ref. As 
an illustration of these results, in Fig. 4 we show differential distributions for Vbb production at 
the LHC for ^/s = 7 TeV, in particular the the of the W boson in Wbb production and the 
invariant mass of the bb pair, in Zbb production. 

These results were obtained in the framework of the Sherpa+OpenLoops Monte Carlo 
generator, using the 5FS, 4FS and doped versions of NNPDF2.3 NLO In each case, a 
consistent treatment of the partonic matrix elements was adopted. From Fig. 4 we see that 
the differential distributions in the doped scheme are typically intermediate between the results 
from the 4F and 5F schemes, as one could expect from the comparison at the level of PDFs 
shown in Figs. 2 and 3. These comparisons imply that the (resummed) change in as(Q) and in 
the PDFs is more significant than the finite-order correction in the matrix element implemented 
in the doped scheme. The results of the calculations in the three schemes are consistent with 
the scale uncertainties of the 5FS calculation. 

Doped PDFs were also used in the ttbb calculation of Ref. where it was found that a 
resummation of the logarithms of PR/rrih arising from the running of as in the doped scheme 
increased to total NLO cross-section by almost 10% as compared to the 4FS calculation. 

The NNPDF3.0 doped PDFs presented in this contribution have been recently used in 
the Sherpa framework to compare the predictions of the 4F, 5F and doped schemes with recent 

















































Figure 4 - Differential distributions for vector boson production with a pair of &-quarks at the LHC 7 TeV: 
the of the W boson in Wbb production (left) and the invariant mass of the bb pair in Zbb production 
(right). The results of the 4F and 5F schemes are compared with the doped scheme. 


ATLAS and CMS measurements at 7 TeV for the differential distributions of Z bosons 
produced in association with 6-jets.^ These comparisons are illustrated in Fig. 5, showing that in 
general the doped scheme leads to a better agreement with the LHC data as compared to the 4F 
and 5F schemes; for instance for the ATLAS yhoost{Z,b) distribution, the doped scheme leads 
to the correct normalization and shape while the 4FS reproduces the shape but undershoots the 
data. 

Summarizing, the doped scheme is a hybrid scheme for the treatment of heavy quark masses 
in perturbative QCD calculations, which has the potential to improve some of the drawbacks of 
4FS calculations, in particular for processes with many colored partons in the final state. The 
NLO and NNLO doped versions of NNPDF3.0 will be made available in LHAPDF6. 
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Figure 5 - Preliminary results for the comparison between the 4FS NLO, doped NLO and 5FS LO 
calculations with Sherpa using the corresponding NNPDF3.0 PDFs, and differential distributions from 
the recent ATLAS and CMS 7 TeV measurements of Z production in association with 6-jets. 
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